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<^ '. ABSTRACT 

O: 

. , Context. The emission of cold dust grains at long wavelengths will soon be observed by the Planck and Herschel satellites and will provide 
, new constraints on the nature of interstellar dust. In particular, the microwave galactic anomalous foreground detected between 10 to 90 GHz, 
proposed as coming from small spinning grains (PAHs), should help to define these species better. Moreover, understanding the fluctuations of 
^ ^ the anomalous foreground quantitatively over the sky is crucial for CMB studies. 

$_( ' Aims. We focus on the long-wavelength emission of interstellar PAHs in their vibrational and rotational transitions. We present here the first 
' model that coherently describes the PAH emission from the near-IR to microwave range. 

I Methods. We take quantum effects into account to describe the rotation of PAHs and compare our results to current models of spinning dust 

■ to assess the validity of the classical treatment used. Between absorptions of stellar photons, we followed the rovibrational radiative cascade of 
PAHs. We used the exact-statistical method of Draine & Li to derive the distribution of PAH internal energy and followed a quantum approach 

^ ■ for the rotational excitation induced by vibrational (IR) transitions. We also examined the influence of the vibrational relaxation scheme and of 
' the low-energy cross-section on the PAH emission. We study the emissivity of spinning PAHs in a variety of physical conditions (radiation field 
, intensity and gas density), search for specific signatures in this emission that can be looked for observationally, and discuss how the anomalous 
T-H ■ foreground may constrain the PAH size distribution. 

CO ' Results. Simultaneously predicting the vibrational and rotational emission of PAHs, our model can explain the observed emission of the Perseus 
molecular cloud from the IR to the microwave range with plausible PAH properties. We show that for A> 3 mm the PAH vibrational emission 
no longer scales with the radiation field intensity (Go), unlike the mid-IR part of the spectrum (which scales with Go). This emission represents 
less than 10% of the total dust emission at 100 GHz. Similarly, we find the broadband emissivity of spinning PAHs per carbon atom to be rather 
constant for Go < 100 and for proton densities < 100 cmr^. In the diffuse ISM, photon exchange and gas-grain interactions play comparable 
roles in exciting the rotation of PAHs, and the emissivity of spinning PAHs is dominated by the contribution of small species (bearing less than 
100 C atoms). We show that the classical description of rotation used in previous works is a good approximation and that unknowns in the 
1 vibrational relaxation scheme and low-energy cross-section affect the PAH rotational emissivity around 30 GHz by less than 15%. 
Conclusions. The contrasted behaviour of the PAH vibrational and rotational emissivities with Go provides a clear prediction that can be 
tested against observations of anomalous and dust mid-IR emissions: this is the subject of a companion paper. Comparison of these emissions 
complemented with radio observations (21 cm or continuum) will provide constraints on the fraction of small species and the electric dipole 
moment of interstellar PAHs. 
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1. Introduction 

The mid-IR spectrum of the interstellar medium (ISM) shows 
prominent bands from 3 to 17 yum, which account for one third 
of the energy emitted by interstellar dust. Such bands are emit- 
ted by very small (subnanometric-sized) dust particles during 
internal energy fluctuations triggered by the absorption of a 
stellar photon (Sellgren 1984). The positions of these bands 
suggest there are arom atic, hydrogenated cycles in th ese grains. 
Leger & Pugetl d 19841) and lAllamandola et aU d 1985k proposed 
polycyclic aromatic hydrocarbons (PAHs) as the carriers of 
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these bands. Despite two decades of experimental and theo- 
retical efforts, the match (band positions and intensities) be- 
tween data only available on small species and observations 
still remains elusiv e, as illustrated recen t ly by the work of 



Peeters et al.l ( l2002h and lKim & Savkallvl ( l2002h . Given their 



important role in the ISM (e.g., gas heating, Habart et al. 2001, 
UV extinction, Joblin et al. 1992), it is necessary to find other 
ways to constrain the properties of interstellar PAHs. 

The Planck and Herschel data will soon trace the emis- 
sion of cold interstellar grains in the interstellar medium. 
Due to their small size, PAHs are heated sporadically (ev- 
ery few months in the diffuse ISM) by absorption of stel- 
lar photons and have a high probability of being in low- 
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energy states (IDraine & Lill200lb . Interstellar PAHs may thus 
contribute significantly to the emission at long wavelength 
(A > 1 mm). In this context, an unexpected emission ex- 
cess called anomalous foreground, correlated t o dust emis 



sion, has been discovered between 1 and 90 GHz jLeitch et al 



1997 ; de OUveira-Costa et al.l2002 ). In this spectral range, sev 
eral galactic emission components (synchrotron, free-free, and 
thermal dust) contribute with a comparable magnitude, and 
only recently has the anomalous foreground been separated in 
WMAP data dMiville-Deschenes et al.ll2008l). Spinning , smal l 
dust grains were first proposed by iDraine & LazarianI (119981) 
(hereafter DL98) as a possible origin to this anomalous com- 
ponent. Since then, analysis of observations has suggested that 
the anomalous foreground is correlat ed to small grain emission 
(lLagachell2003 : ICasassus et al ] l2006h . 

In this paper, we study the emission of PAHs with particu- 
lar emphasis on the long-wavelength part of the spectrum. In 
this spectral range, the emission is dominated by species in 
low-energy states for which each photon exchange represents a 
large energy fluctuation. To derive the internal energy distribu- 
ti on of PAHs, we use h ere the exact statistical method described 
in Draine (feLil ( 2001 ). We include low-frequency bands to the 
vibrational mode spectrum of interstellar PAHs and examine 
the influence of the internal vibrational redistribution hypothe- 
sis (see Sect. [3]) often used to describe their vibrational relax- 
ation. We describe the rotation of PAHs with a quantum ap- 
proach where specific processes in the rovibrational relaxation 
are naturally included (molecular recoil after photon emission; 
rovibrational transitions that leave the angular momentum of 
the molecule unchanged or Q bands). We present in a variety 
of interstellar phases the rovibrational (IR) and broadband rota- 
tional (spinning) emission of PAHs in a consistent fashion. Our 
model results are compared to observations of the interstellar 
emission from the IR to the microwave range and show how 
the size distribution and electric dipole moments of interstellar 
PAHs may be constrained. 

The paper is organized as follows. Section 2 describes the prop- 
erties of PAHs adopted in this work. Section 3 discusses the in- 
ternal energy distribution of PAHs and associated rovibrational 
emission. Sections 4 and 5 discuss the rotational excitation and 
emission of PAHs. In Sect. 6, we apply our model to the case 
of a molecular cloud in the Perseus arm. Finally, conclusions 
and observational perspectives are given in Sect. 7. 

2. The properties of interstellar PAHs 

In the ISM, the stablest PAHs are found to be compact species 
dLeger et alj|l989l: iLe Page et al.ll2003l) . SmaU PAHs are pla- 
nar, but above some poorly known size threshold (40 to 
100 C atoms), interstellar formation routes may favour three- 
dimensional species (bowl- or cage-shaped) containing pentag- 
onal cycles (see Moutou et al. 2000 and references therein). As 
we see later, small PAHs (containing less than 100 C atoms) 
dominate the rotational emission, so we assume that interstel- 
lar PAHs are planar with a hexagonal (D^f,) symmetry. The 
PAH radius is thus a(A) = 0.9 VA^, where N c is the num- 
ber of carbon atoms in the grain (IOmonllll986l) . The formula 
for such molecules is C^piH^p, and their hydrogen-to-carbon 



ratio is ^ - /h-^J^, where /h is the hydrogenation fraction 
of PAHs. In this work, we assume that PAHs are fully hydro- 
genated (/h = 1). 



2.1. Absorption cross-section 

The excitation, cooling, and emission of PAHs depend on their 
absorption cross-section (Sect. |3]l, which we describe now. 
We too k the visible-UV cross-section from lVerstraete & Legei 
(119921) and applied their size-dependent cut-off for electronic 
transitions in the visible-NIR range. Th e resulting cro ss-section 



compares well to the available data dJoblin et alj Il992i) . In 



Appendix |A] we discuss the mid-IR bands considered in this 
work. Each vibrational mode is assumed to be h armonic, and 
the co rresponding band profile has a Drude shape ( Draine & Lil 



2001h . The width Av is inferred from astronomical spectra, and 



the peak value cr is chosen so that the integrated cross-section 
crAv is equal to the value measured in the laboratory. By adopt- 
ing the observed bandwidth, we empirically acco unt for the 



complex molecular relaxation and band broadening (IPech et al 



12002; Mulas_et aL 2006) in i ntersteUar PAHs. We adopt here 



the band strengths given in iPech et al.l (l2002h that were in 



ferred from laboratory data. Different d efinitions of the PAH IR 
cross-sections have been proposed by iRapacioli et al.l (120051) . 
Flagev et alJ (l2006l) . and Draine & Li (2007). 



We used the database of lMalloci et al.l ( 120071) to define an aver- 
age broadband cross-section of the far-IR vibrations of PAHs. 
At frequencies below 500 cm"', each species features many 
modes. However, for compact species, modes accumulate in 
three definite frequency ranges: modes with a frequency of 
less than 100 cm"', modes between 100 and 200 cm"', and 
modes between 200 and 500 cm"'. We therefore model the 
far-IR cross-section of compact PAHs with 3 modes (Table 
[T]). The frequency of each mode is the average of all modes 
falling within the given energy range, weighted by their corre- 
sponding integrated cross-sections. We find that the frequency 
of the lowest energy mode depends on the molecular size as 
A^^;' (see Fig. [I])- Conversely, for the two other modes, the av- 
erage energy is instead independent of the size. The integrated 
cross-section for these 3 modes was estimated as follows. From 
the Mallocci database, we first derived the fraction of crAv for 
each of the 3 modes. Then, we assumed that the total crAv be- 
low 500 c m"' is given b y the in tegral of the absorption cross- 
section of ISchutte et"aD (Il993h : ctfir = 4.3 x 10"^" ^ ' "^ 
cm^ per C-atom. The integrated cross-sections of each of the 
3 modes was finally obtained by multiplying the former value 
by the crAv fractions inferred from the database. 

The parameters of the far-IR bands adopted here are given in 
Table [T] and the full cross-section is displayed in Fig. |2l We 
show in Sect.|3]that this set of IR bands provides a good match 
to observed interstellar spectra. Finally, we note that the in- 
plane or out-of -plane character of each vibrational band is im- 
portant because of the different associated weights in the ro- 
tational excitation (Sect. 14.1b . This character is indicated in 



Tables lATTI for the mid-IR bands dSocratesluOOll) . In the case 
of the far-IR bands, this character is not as well known and 
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Table 1. Far-IR rovibrational bands of PAH cations adopted in 
this work, with the percentage of the total oscillator strength in 
each band. 
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° Similar bands were obtained for PAH neutrals. 

* In-plane (ip) or out-of-plane (op) character of the band. 
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Fig. 1. Frequencies of the first (lowest energy) vibrational mode 
versus Nc for PAH cations (black diamonds) and PAH neutrals 
(grey circles) from the Mallocci databas^see text). The solid 
line is the relationship we adopt between the band position of 
cations and Nc, and the dashed line shows the case of neutrals. 
The dotted lines show extreme cases for this relationship. 



Fig. 2. Absorption cross-section of PAH cations per carbon 
atom (1 Mb = IQ-"^ cm^). The solid, dotted, dashed, and dot- 
dashed lines show the cases for Nc = 24, 54, 96, and 216, re- 
spectively. 
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Fig. 3. Symmetric top molecule: Da, Db, and Dc are the prin- 
cipal axis of inertia and J the total angular momentum of the 
molecule with K its projection along Dc- 



we assume that 1/3 (2/3) of the oscillator strength come from 
out-of-plane (in-plane) transitions respectively (Table [l]!. 

2.2. The rigid rotor model 

While describing the rotation of a molecule, the relevant oper- 
ator is the total angular momentum 7, which includes the elec- 
trons and nuclei contributions without the spin. We note that 
Da, Db, and Dc are the principal axis of inertia. We assumed 
that PAHs are oblate symmetric top molecules with the axis Dc 
perpendicular to the plane of the molecule and parallel to Oz- 
We called Ic the inertia moment with respect to Dc and Ia, Ib 
the inertia moments with respect to Da and Db, which were 



- The molecules considered are perylene {CiqHh), 
benzo[g,h,i]perylene (C21H12), coronene (C24-f/i2), bisanthene 
(C28^^i4), dibenzo[bc,ef]coronene (Cjo-Wm), ovalene (C32-f/i4), cir- 
cumbiphenyl (Csg/fis), circumanthracene (C^oHif,), cirumpyrene 
(C42//16), and circumovalene (C(,f,H2o). 



taken to be parallel to Ox and Oy. The rotational Hamiltonian 
is then 

Jl A' 

H = -i- -H — -H — (1) 

2Ia 2Ib 2Ic 

where 7,, 7,, and J, are the projections of ~I along the three 
inertia axis. Given the large number of carbon atoms in inter- 
stellar PAHs, we assume here that they are uniform disks with 
Ia - Ib - = (symmetric tops), where M is the molecu- 
lar mass. With Dc \\ Oz, we get 

and the rotational energy is 

Ero,^BJ(J+l) + (C-B)K^ (3) 

where B/hc = - 7 A^^^- cm-' (neglecting the contribu- 

tions of H to the molecular mass) and C = B/2 are the rota- 
tional constants associated to Db and Dc- The quantum number 
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K is the absolute value of the /--eigenvalues. For a symmet- 
ric top molecule, the selection rules for rovibrational electric 
dipole transitions are A/ = 0, +1 and /S.K - 0, +1. By defining 
the dipole moment as 7? = J? + y Ug, with jt, along Dc and ju g 
in the molecular plane along Dg ( Townes & Schawlowlll975 ). 
two kinds of transitions can be distinguished. First, the parallel 
transitions with IS.K - for which the change of dipole moment 
in the transition is parallel to the top axis {Dc) of the molecule. 
Second, the perpendicular transitions with AK = ±1, for which 
the change of dipole moment is perpendicular to Dc- Parallel 
transitions thus correspond to out-of-plane vibrations, whereas 
perpendicular transitions correspond to in-plane vibrations. 
Since the available microwave data are broadband observations 
(/i/A/l of the order of a few), we make several simplifying as- 
sumptions in the description of the rotational motion of PAHs. 
We thus assume that the rotational constant B is the same in 
all vibrational levels, and within the framework of a rigid rotor 
model, we neglect the centrifugal distortion terms in t he energy 
equation that are usually s mall for large molecules ( Herzberg 
1968abl:lLovas et al.ll2005l) . 



2.3. Electric dipole moment 

The rotational emission of PAHs depends on their permanent 
electric dipole moment, /i. Symmetric (De/,), neutral, and fully 
hydrogenated PAHs have /i ~ 0. Spectroscopic analysis of their 
IR emission bands suggests that interstellar P AHs can hqve 



a cati onic form that is partially h ydrogenated jLe Page et af. 
2003h . and maybe also substituted jPeeters et al.ll2002l l2004ar 
For instance, a PAH having lost one H atom has p ~ 0.8 to 1 
D, and a PAH cation where a C atom has been substituted by 
N would also have // ~ 0.1 to 1.5 D depending on its size (T. 
Pino, private communication). Moreover, it has been proposed 
that non-planar PAHs containing pentagonal rings may exist in 
the ISM (see Moutou et al. 2000 and references therein): such 
species are known to have large dipole m oments as re cently 
measured on coranulene, €20^/^10, p - 2D jLovas et al ]|2005). 
In this work, we express the electric dipole moment of inter- 
stellar PAHs as in DL98: 

p(Nc) = m y/K, + 4.3 X 10"^ Va^Z ^ m (Debye) (4) 

where Nai is the total atoms number in the molecule, Z is its 
charge, and m a constant. Unless otherwise stated, we will use 
m = 0.4 D. 



3. Internal energy distribution and rovibrational IR 
emission of isolated interstellar PAHs 

After the absorption of a visible-UV photon, a PAH cools off 
by emitting IR rovibrational photons. These photons reduce its 
angular momentum and also may increase it by recoil (a purely 
quantum effect, see Sect. 14. 11 1. Previous studies have mostly 
used a thermal description of molecular cooling. In fact, since 
PAHs spend a large fraction of their time at low internal ener- 
gies (see Fig.|4]i, their emission at long wavelengths cannot be 
considered as a negligible fluctuation compared to their total 
internal energy. While estimating the emission of PAHs at low 
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Fig. 4. Internal energy distribution of PAHs with Nc - 24 
(solid line and box for P(0)) and Nc =216 (dashed line and tri- 
angle) in the case of a) an MC (Go = 10"^), b) the diff'use ISM 
(Go = 1) and c) the Orion Bar (the radiation field is the sum of 
the CMB, the ISRF and a blackbody at 37 000 K corresponding 
to Go = 14, 000). To illustrate the effect of the radiation field 
hardness, we show in b) the case of an Orion Bar type radiation 
field scaled down to Go = 1 (grey lines). 



temperatures, the validity of a thermal approach is question- 
able. Because of the rapid energy redistribution between inter- 
actions (photon exchange or collision with gas phase species), 
PAHs rapidly reach the rmodynamical equilibrium while iso- 



lated (ILeger et al.llI989h . To describe this situat ion, we used 



the exact-statistical method of iDraine & Lil (1200 1 h to derive the 
stationary internal energy distribution of PAHs. 

The energy distribution, P{E), depends on (a) the energy den- 
sity of the exciting radiation field ug - A-nvIy (where is the 
brightness), (b) the absorption cross-section, cr, of interstellar 
PAHs, and (c) their rovibrational density of states, p(E). For the 
stellar contribution, we used the local interstellar radiation field 
(ISRF, Mathis et al. 1983) or a blackbody. We sometimes scale 
this stellar radiation field with a factor. Go, equal to I in the case 
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of theMathis fieldd Throughout this work, the contribution of 
the CMB is included in the radiation field. The vibrational ab- 
sorption cross-section used has been described previously, and 
for simplicity, we did not include rotational bands. Unless oth- 
erwise stated, we took the cross-section of PAH cations. The 
density of states was obtained by first deriving the vibrational 
mode spec trum from a Debye model and then applying the al- 
gorithm of Beyer & Swinehard (Il973h for each molecular size 
(see Appendix iBt. 

After each absorption of a stellar photon, we assume that the 
excitation energy of the PAH is rapidly re distributed among 
all the vibrational modes ( Mulas et alJ2006l) . This radiationless 
and isoenergetic process is called internal vibrational redistri- 
bution (IVR). The energy distribution is then com puted accord 



200 Ih 



ing to the exact-statistical formalism described in iDraine & Li 



( 1200 11) . The energy scale (referenced to the zero point energy) 
is divided into bins of energy £, and width AE,- with / = 
0, 1, ...,M (M = 500) with < < 2El where = 1.110^ 
cm'' is the Lyman limit. When its energy is below that of the 
first excited vibrational state, the molecule is in a rotational 
state. Letting P, be the probability of having a PAH in the en- 
ergy bin / we have 



dPi 

It ~ /-i 



(5) 



where Tj^i is the transition rate from the state j to the state / for 
a molecule. We solve thi s equation in t he stat ionary case. More 
details can be found in Draine &"l3 (2001). Figure H] shows 
the energy distribution for three different radiation fields, repre- 
sentative of a molecular cloud, the diffuse interstellar medium 
(DIM), and the Orion Bar We see that the most probable en- 
ergy increases when the radiation field intensity Go increases, 
as expected for PAHs whose cooling by IR emission is inter- 
rupted more frequently by absorption events. The sharp cut-ofF 
at ~ 1.1 X lO'' cm"' stems from the Lyman limit of photon 
energies in neutral regions. We also observe a tail at higher 
energies that comes from to multiphoton events, i.e., absorp- 
tion of a second photon while the PAH has not completely 
cooled off. This tail becomes more significant as absorption 
events are more frequent, i.e., for large PAHs or intense radi- 
ation fields. Conversely, we show in Fig. that the radiation 
field hardness has little influence on P(E). To estimate the ef- 
fect of the rotational absorption that we neglected, we included 
a band centred at 1 cm ' of width 1 cm ' and corresponding to 
J = 150 (a-INc ~ 9 X 10^22 cm^/C, see Sects. |42]and|43]i: we 
found P{E) to be twice as large between 1 and 10 cm ' and un- 
changed otherwise. This would change the emission around 1 
cm at very low flux levels. The rotational excitation rates (Sect. 
I4.2l i. which depend on the populations of excited vibrational 
levels, are unaffected by this hypothesis. 
Knowing the internal energy distribution of PAHs, we can de- 
duce their IR emission from the upper state u by summing 
the contributions of all lower vibrational modes / (IDraine & Li 



^ Go scales the radiation field intensity integrated between 6 and 
13.6 eV. The Mathis radiation field, Go = 1, corresponds to an inten- 
sity of L6 X 10"^ erg/s/cm^. 



2hv I uf \ 

.^,._.(,)P(,,)(l,_£_) 



(6) 
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where P{hv) = V P, V ^ A£„ G„/(/iv) 

^ g, 



I 11=0 



(7) 



is the number of rovibrational p hotons emitted at e nergy liv. 
The Gui functions are defined in IDraine & Lil ( 1200 lb . The de- 
generacies g„ and gi are the numbers of energy states in bins u 
and I, respectively (see AppendixlBj. Figure|5]shows a compar- 
ison between mid-IR observations and our mode l results with 
a pow er law size distribution ns(a) ~ a'^'^ (as in iMathis et al. 
(Il977h hereafter M RN) that provid es a simple representation 
of the actual «j(a) (IKim et al.lll994i) . We see that the IR emis- 



sion from the DIM is explained well with a standard abun- 
d ance of carbon in P AHs and as large fraction of cations as 
in lFlagev et al. I (l2006l) . A similarly good match is obtained for 



the Orion Bar spectrum with strong PAH depletion probably 
reflecting an efficient photodestruction in this excited environ- 
ment. Figure [6^ shows the behaviour of the PAH vibrational 
emission with the radiation field intensity Go: whereas the IR 
part scales with Go, the mm-range {A > 3 mm) does not. The 
emission in this spectral range is produced mostly by PAHs 
in low-energy states for which the dominant heating source is 
absorption of CMB photons, and this is why this part of the 
spectrum does not vary with Go. Finally, we emphasize that 
at low energies {A > 20 fim), our model shows the broad- 
band behaviour of the PAH vibrational emission. When a de- 
tailed vibrational mode spectrum is used, the PAH emission 
at /I > 20 jum is a superposition of numerous narrow bands 
that may be detectable with some of the instruments onboard 
Herschel dMulas et al.ll2006l) . 



3.1. Decoupling of vibrational modes 

Our derivation of P{E) assumes an efficient energy redistribu- 
tion between vibrational modes (IVR) during the PAH relax- 
ation. The IVR thus involves a coupling between vibrational 
modes via intramolecular transitions. However, it is known 
that IVR is no longer efficient when the e xcitation energy o f 
the molecule is below some threshold E^ec dMulas et al ]l2006l) . 
and we will then speak of IVR breakdown. When E < E^ec, 
the excitation of energetically accessible vibrational modes is 
frozen according to microcanonical statistics at energy E^ec- 
The molecule then continues to cool via emission of the excited 
vibrational modes. To study the influence of decoupling on the 
long-wavelength emission, we consider two extreme cases for 
the relaxation at £ < E^ec'- 

- case 1: no IVR, cooling by IR forbidden modes. We assume 
these modes to be defined like those of Table [T] but with an 
oscillator strength 10"* lower and, 

- case 2: no IVR, cooling only by the first vibrational mode at 
2040/A^c (cm-i) as defined in Table [T] 
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Fig. 5. Left: Dust IR emission of the DIM (solid black lines and symbols) with A^// ~ 10 cm ( lBoulangej2000h . The model dust 



emission is overlaid in bold grey. MRN type size distributions are used for all dust populations. The R^H vibrational contribution 
(dashed line) is from the present model with Nq - [18,96], a 1:3 mixture of neutrals and cations, full hydrog en coverage, and 
43 ppm of carbon. The contributions of larger grains are from the model described in lCompiegne et alj ( 120081) . The dotted line 
is the contribution of graphitic very small grains of radii a - 0.9 to 4 nm containing 39 ppm of C. The dot-dashed line is the 
contribution of silicate big grains with a - 0.4 to 250 nm using 37 ppm of Si. Right: Dust emission from the Orion Bar (the noisy 
black line is the ISO-SWS spectrum) with Nh ~ 1.8 x 10^' cm"^. The model dust parameters are as above but for PAHs less 
abundant (18 ppm of C) and more ionized (82% are cations). 
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Fig. 6. PAH rovibrational emissivity for a MRN size distribu- 
tion with Nc = 18 - 216. In a) we show the case of radiation 
field intensities corresponding to Go = 1 and 100 and in b) we 
illustrate the effect of IVR breakdown in case 1 (dotted line) 
and case 2 (dot-dashed line) (see section lTTT i. 



and compare them to our model with no decoupling (IVR al- 
ways fulfilled). We take Ejec to be constant and equal to 500 
cm"', a value relevant for P AHs containing 20 to 30 carbon 
atoms ( Joblin & Mulasll2008 ). Actually, E^ec is expected to de- 
crease with molecular size b ecause of the increasing density 
of states at a given energy ( IMulasI 119981) : the former cases 
therefore provide an upper limit to the variations in the PAH 
long-wavelength emissivity. FigurejSJ? shows that, as expected, 
case 1 provides a lower limit to the long-wavelength emis- 
sion of PAHs, whereas case 2 is an upper Umit. The model 



assuming IVR is close to case 2. A comparison of our model 
band fluxes (wi t h IVR ) with the Monte-Carlo simulations of 
Joblin & MulasI (l2008h . in the case of the coronene cation 
(C24HI2), shows good general agreement, especially for the 
low-energy modes. We conclude that the IVR hypothesis pro- 
vides an upper limit to the PAH vibrational emission at A > 3 
mm and that the conclusions of the former section hold. 



4. Angular momentum distribution 

Building the angular momentum distribution of interstellar 
PAHs is driven by photon exchanges and gas-grain interactions 
(DL98). In the DIM, pervaded by the ISRF, and for a PAH 
bearing 50 carbon atoms, the mean time between absorptions 
of visible-UV photons (~ 0.2 yr) is comparable to the mean 
time between emissions of rotational photons, as well as to the 
mean time between PAH-hydrogen collisions (for a gas density 
of 100 cm"^ and a temperature of 100 K). Photons absorbed 
in the visible-UV have a weak eff'ect on the total angular mo- 
mentum; indeed, each photon exchanged carries a unit angular 
momentum, and the numerou^ IR photons emitted overwhelm 
the influence of the photon absorbed. To estimate the angular 
momentum distribution, we take the following processes into 
account: 

- IR rovibrational photon emission 

- purely rotational photon emission 

- H2 formation 

- collisions with gas (neutral atoms and ions) 

- plasma drag 



* Energy conservation implies tliat 40 IR phiotons are emitted after 
each absorption. 
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- photoelectric effect 

All these processes lead to a change in the angular momentum 
of the molecule, some of them excite the rotation, and others 
damp it. For the radiative processes we adopt a quantum ap- 
proach where each rotational state is tre ated individually an d 
transition rules applied. As discussed by iRouan et al.l (119971) . 
the building of the angular momentum distribution, n(J), can 
be considered as a stationary random walk in a potential well 
with a minimum for J - Jq, for which the rate of /-change is 
zero: 







(8) 



'Jo 



with AJ the change of / produced by the eve nt number 
and T the mean time between two events /. As in Rouan et al 



(Il992h . we assume an efficient intramolecular vibration to rota- 
tion energy transfer (IVRET). We take «(/) to be the same for 
all vibrational levels, equal to a Maxwell distribution n{J) = 
noJ^ exp(-J^ / J^), where »n is a norm a lizatio n factor. Indeed, 
MulasI d 19981) and lAli-Haimoud et al. (l2009l) show that this 
form of n(J) is a good approximatior 1 In the following, we 
establish the rate of /-change due to the rovibrational emission 
of isolated PAHs modelled in Sect. [3] We also present the rate 
of purely rotational photon emission, while the contribution of 
the gas-grain interactions is described in AppendixlP] 

4.1. Rovibrational transitions 

The rovibrational IR emission can be both an exciting and a 
damping process for the rotation of a PAH. Assuming that 
the interstellar PAHs are symmetric top molecules, the selec- 
tion rules for the emission of an IR rovibrational photon are: 
AJ = 0, +1, and AK = 0, +1. Transitions corresponding to 
AJ - +1,-1, and are called the P, Q, and /?-bands, re- 
spectively. As seen in Sect. 12.21 two types of transitions are 
possible: parallel ones (out-of-plane vibrational motion) with 
AK = and perpendicular ones (in-plane vibrational mo- 
tion) with AK - +1. The type of each transition is given in 
Tables lA.ll and [T] For a rovibrational transition, (v, J, K) — > 
(v - 1,7 + AJ, K + AK), the rate is proportional to the sponta- 
neous emission coefficient and to the probability for the grain 
to lose the corresponding transition energy. The spontaneous 
emission rate is proportional to the Akj factors, which repre- 
sent the angular part of the transition probability. Formulae for 
these factors are given in Appendix |C] Expressed in terms of 
cross-section, the transition rate is (with v, and Av, in cm '): 

Wfl' = Sttc (cr,Av,) i^ff ^l^y P^hv^") (9) 

where the Wf^^'' give the transition rate for AJ - 0, ±1, and 
for the vibrational mode ; at a frequency vj^'". This frequency 
depends on v,o (the frequency of the vibrational transition for 
J - Q) and on J and K (see Appendix O- Since the /- and 



^ With this form of n{J) it is possible to define a rotational temper- 
ature (see Appendix lEt. 



/T-terms, in v^^'^, are always much smaller than v,o, we take 
P(hvf^'^) - P(hvio). We note that the in-plane transitions pro- 
vide higher (by a factor ~ 2) rates than the out-of-plane ones, as 
a consequence of the Akj sum values. The total transition rates 
are obtained by summing over all bands: = 2, Wf^'^. 

Finally, the rate of change in / due to rovibrational transitions 
is (T"'Ay)/R = - W'. Rotational excitation dominates at 
low /, whereas damping is the dominant process for high J- 
values. For vibrational modes at higher frequency, - 
is smaller because the rotational terms are smaller in v,- (see 
App.O: (t 'A7)/r is therefore dominated by the contribution 
of vibrational modes at low-frequencies. 

4.2. Rotational emission 

We consider here the spontaneous emission of purely rotational 
photons. The selection rules for such tra nsitions are A/ = -1 
and AK = jTownes & Schawlowlll975h . No change in K oc- 
curs because, for rotational transitions, the dipole moment of 
a symmetric top molecule necessarily lies along its symmetry 
axis. The transition rate is then simply related to the sponta- 
neous emission coefficient Ay y_ i : 



(t-'aju, 

Ajj-i 



-Aj,j-i 
5Uk^ 



bVj' 



(2J +\)- -{J + 2) 
3(27+1)2 



(10) 
(11) 



where fi is the electric dipole moment of the molecule and the 
factor Akj- has been used (Appendix 0. In the high-/ limit, 
hvxA J J -I tends to the classical expression of Larmor. Finally, 
the damping rate due to rotational emission i^ 



{t-'AJ\„, 



^3 (27 +1)2 -(7 + 2) _i 

X J s 

(27+ 1)2 



(12) 



4.3. Equilibrium angular momentum Jq 

We estimate damping and excitation rotational rates for the in- 
terstellar phases described in Table |2] When solving Eq.[8] we 
obtain the equilibrium angular momentum 7o. Figure [T] shows 
7o for the DIM and the Orion Bar. In all cases, 7o increases with 
Nc because there are more steps in the random walk; indeed, 
the cross-section and absorption rate scale with Nc (see Fig. 
|2|, so more IR photons are emitted for larger A^^. Similarly, 7o 
reaches higher values when the radiation field intensity or the 
gas density increases (see Fig.[8]l. 

In the DIM, the excitation comes from gas-grain collisions, 
whereas the damping is dominated by the emission of rota- 
tional photons for Nc < 140 and by IR photon emission for 
larger PAHs ((t"'A7),o/ N^f', whereas (t"'A7)/r oc Nc)- We 
note similar 7o-values for the CNM, WNM and WIM (Fig.|7h) 
in spite of the large differences in gas density and temperature 



^ In the limit 7 » 1 , we estimated that the stimulated emission and 
absorption of CMB photons represent less than 20% of the sponta- 
neous emission rate. 
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Table 2. Physical parameters for the typical interstellar phases 
considered, nu the hydrogen density, T the gas temperature, 
the electrons density, «//+ the proton density, and nc+ the density 
of atomic carbon ions. 



(nm) 







CNM^ 


WNIVF 


WIM^ 


Orion Bar 


Go 


0.01 


1 


1 


1 


14 000 




300 


30 


0.4 


0.1 


10^ 


r(K) 


20 


100 


6000 


8000 


400 


He (cm^^) 


0.03 


0.045 


0.04 


0.1 


3 







1.2 X 10-3 


0.1 


0.99 


10-^ 


"C+ /"H 


lo-"* 


1.3 X 10-"* 


1.3 X 10^" 


10-3 


2 X lO-'* 



" Molecular cloud. 
* Cold neutral medium. 
' Warm neutral medium. 
Warm ionized medium. 



(see Table |2]i. In the Orion Bar, because of the intense radia- 
tion field, the excitation is driven by the photoelectric effect, 
whereas the damping is dominated by the emission of IR pho- 
tons. For the case of MCs (not shown), Jq also rises and reaches 
high values (/q - 1300 for Nq - 200). Indeed, the damping by 
IR photon emission is no longer efficient (Go ~ 0.01), while the 
gas-grain rates are strong because of higher density. We show 
in Fig. |7j3 the influence of different choices of the vibrational 
mode spectrum (or cross-section) on /o: changes are below 
20% and affect only the large sizes {Nq > 100), which make a 
minor contribution to the rotational emissivity (see Sect.|5]i. We 
found comparable variations of Jq for the vibrational relaxation 
(IVR) cases 1 and 2, discussed in Sect. 13. II 
The processes contributing to the rotation of PAHs depend on 
the intensity of the radiation field (Go) and on the gas density 
(nn). We now examine the influence of these parameters on Jq. 
Other important quantities (gas temperature; electron, proton, 
and ab undances) are obtain ed at thermal equilibrium with 
CLOUDY ( lFerland etalJll998H1 . To study the influence of Go, 
we took the radiation field to be a blackbody of effective tem- 
perature 22,000 K and scaled it in order to vary Go between 
10^^ and 10^. The shape of the radiation field is then always 
the sam^^. Figure [8^ shows Jq as a function of Go for different 
PAH sizes. As discussed before, Jq increases with Go because 
the excitation rate from IR emission (which scales with Go) is 
dominant. This rise becomes steeper for larger PAHs because 
the frequency of the first vibrational mode decreases (as A^^') 
and requires lower values of Go to be excited (the maximum 
of P(E) moves to higher energies as Go increases, see Fig.H]:). 
However, we note that Jo is affected little by variations of Go 
over the range 0.01-100. The result of varying hh is shown in 
Fig. [Sji. In all cases, the incident radiation field is a blackbody 
with Teff - 22,000 K and Go = 1 . The gas-grain processes be- 
come dominant for «// > 30 cm"-' and Nc > 50. At lower gas 



' We assumed 130 (320) ppm of carbon (oxygen) to be in the gas 
phase. Parameters were taken from the optically thin zone of isochoric 
simulations with CLOUDY. 

We also varied the shape of the blackbody keeping Go constant 
(as well as nu). Similar Jq values were found for r^// between 10* 
and 5 x 10* K. 



l.OxlO' 




Fig. 7. Jo versus Nc for PAH cations. In a) we show the case of 
the diffuse medium. Solid line shows the CNM, dashed line the 
WNM and dotted line the WIM. In b) we illustrate the influence 
of the absorption cross-section in the CNM: cations (solid line), 
neutrals (dashed line) and extreme behaviours (grey lines) for 
the first vibrational mode versus Nc (see Fig. [Til are shown. In 
c), we present the case of the Orion Bar. 



densities, radiative processes prevail (excitation by IR emission 
and damping by rotational emission), and Jo does not change. 
We therefore expect the rotational excitat ion of small PAHs t o 
hardly be variable in the DIM as found bv lDavies et al. I (l2006h . 



4.4. Modelling the rotational excitation by IR emission 

Following the absorption of a stellar visible-UV photon, the 
rovibrational cascade of PAHs is a complex process that in- 
volves many states. Molecular, state-to-state models require a 
detailed database, use Monte-Carlo simulations (see Mulas et 
al. 2006 for references), and so far do not include the gas-grain 
interactions that are important. Given our incomplete knowl- 
edge of interstellar PAHs and the fact that their rotational emis- 
sion may so far be seen in broadband data, former models of 
spinning dust (DL98 and Ali-Haimoud et al. 2009) made sim- 
plifying assumptions to describe radiative processes and per- 
formed a classical treatment of the gas-grain interactions]^. We 
assess here the impact of these assumptions on the rotational 
emissivity of PAHs. 

First, the internal energy distribution of PAHs was derived in 
the thermal approximation. As discussed in Sect. IH this is a 
questionable assumption for describing the long-wavelength 



' These models can also be applied to describe the rotational emis- 
sion of other grain types. 
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Fig. 8. a) Effect of varying Go on Jo for PAH cations with Nc = 20, 30, 40, and 100 (from bottom to top). The values for A^c = 100 
(triangles) have been divided by 2. The gas density is 30 cm"^ and the other physical parameters for the gas are determined using 
CLOUDY ( Ferland et al.ll 1998 ). b) Same as a) but for varying «// at Gq - I. In both cases the radiation field is a blackbody with 
Teff = 22 000 K and the gas parameters have been obtained with CLOUDY at thermal equilibrium (see text). 
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Fig. 9. Rate of /-change due to rovibrational emission for a 
PAH with Nc - 24 and heated by the ISRF. a) Comparison 
with former models: o ur rate (so l id line ), the rate of DL98 cor- 



rected as in lAli-Hai'moud et al 



rate of lAU-Haimoud et al.l(l2009l) (dotted line), h) Effect of IVR 



(l2009l) (dashed line) and the 



breakdown, case 1 and 2 (see Sect. 13. 11 1. 



emission of PAHs and the change in angular momentum it in- 
duces. Moreover, the excitation rate by IR emission (the re- 
coil due to emission of individual photons), which is a purely 
quantum effect, has been described as a random walk of the an- 
gular momentum starting with a non-rotating grair0- Finally, 
an efficient vibrational redistribution (IVR, see Sect. 13.11 ) was 
assumed throughout the energy cascade following a photon ab- 
sorption. In this work, we improve on the first aspect by de- 
riving the internal energy distribution of isolated PAHs using a 
microcanonical formalism and including the rotational density 
of states. Next, we follow a quantum approach to treat the rovi- 



Conversely, the damping part is correctly described by a classical 
model when 7 >> 1 (Ah-Haimoud et al. 2009). 



brational emission where the recoil due to photon emission and 
selection rules are naturally included. In Fig. |9^, we compare 
the absolute value of our rate of angular momentum change by 
rovibrational IR emission |(T"'Ay)/R| - - (see Sect. 
14. It to former works. For a PAH containing 24 carbon atoms, 
all rates decrease with the angular momentum and cross zero 
for J between 100 and 150 (singular points in our logarithmic 
representation). The main difference with previous models is 
that the excitation-to-damping transition (the zero value) oc- 
curs at higher /-values in our case. This discrepancy dimin- 
ishes with increasing sizes, and identical IR rates are found for 
species with more than 100 C atoms. 

We show in Fig. |9j' the effect of relaxation schemes other than 
IVR. Between extreme cases 1 and 2 (see Sect. 13. lb . the damp- 
ing IR rate is multiplied by a factor 2, which leads to a bal- 
ance at a lower /-value. Variations in (t"' Ay)/R are in fact quite 
similar to those induced by the other approximations discussed 
above. All these variations affect the rotational emissivity of 
PAHs, around 30 GHz, by at most 15% (in the CNM and in- 
cluding the gas-grain interactions^ We conclude that, in spite 
of the assumptions made, former spinning dust models (which 
are fast computationally) provide a sufficiently accurate rota- 
tional emissivity. 

5. Spinning dust emission 

The power emitted by a PAH containing Nq carbon atoms in a 
rotational transition from state / to / - 1 is equal to 



(13) 



when taking the spontaneous emission rate A y y_i from Eq.fTOl 
and with the transition energy hv = 2BJ (AK = for rota- 
tional transitions). With B oc A^^^ and fr oc A^^, we find that 



" A correction of the same order was found by lAli-Hai'moud et all 
i l2009l) . who used the Fokker-Planck equation to derive n(J) instead of 
a Maxwell distribution. 
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Fig. 10. Rotational emissivity of PAH cations. We take 430 ppm of carbon in PAHs, Nc =24-216 and assume m - 0.4 Debye. 
Panel a): effect of changing the fraction of large PAHs (decreasing with a = to 6). The case of a log-normal distribution 
centered around Nc - 44 and width 0.4 is also shown. Panel b): rotational emission spectrum in different environments with a 
power law size distribution and a = 3.5. 



P(J) oc A^^^: small PAHs will therefore have a dominant con- 
tribution to the rotational emission. For instance, the peak value 
of the rotational emissivity of a PAH with Nc = 96 is 10 times 
lower than for a species with Nc - 24 (see also Ali-Haimoud 
et al. 2009). With n ,, the size distribution of PAHs (the number 
of PAHs of a given size a or Nc per proton), and the angular 
momentum distribution n{Nc, J), we get 



J y - 

4;r jM 



' Ay^y_, n(Nc,J) ^ nANc) dNc 
2Bc 



(14) 



Since we are interested in a broadband spectrum, we take the 
rotational bandwidth to be 2Bc. As discussed before, we as- 
sume that n{Nc, J) is described correctly by a Maxwell dis- 
tribution. Given the relationship Jo{Nc) (Fig. |7]i and the size 
distribution, the emission of spinning PAHs can be calculated 
from Eq.[T4l It scales with x S pah- We illustrate the rela- 
tive significance of small species in Fig. [TOh by changing the 
index of the power law size distribution, n.s{a) ~ a'". Defining 
small PAHs by Nc < 100, they represent 50% (97%) of the to- 
tal abundance for a - (6). As expected, the rotational emis- 
sivity for ff = 6 is highest. It is also blueshifted by 5 GHz 
and broadened by about the same amount with respect to the 
a = case. Figure [TOb shows the case of different interstellar 
environments and illustrate the influence of the gas-grain pro- 
cesses, which become strong in the CNM (while subdominant 
in the WIM/WNM). 

We emphasize here that the level, peak position and width 
of the rotational emission spectrum depend on (a) the frac- 
tion of small PAHs {Nc < 100), (b) m the scaling factor for 
the dipole moment, and (c) the physical parameters of the gas 
along the line of sight (temperature, density, and ionization). 
Microwave observations alone cannot constrain all these un- 
knowns. A quantitative description of the spinning dust emis- 
sion will require near-IR data (to constrain the fraction of small 
PAHs) and radio maps (21 cm and continuum) to derive the 
physical state of the gas. 
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Fig. 11. Comparison to observations (symbols) of the Perseus 
molecular cloud G159.6-18.5. The dashed lines show the free- 
free (low frequency) and big grains (high-frequency) emis- 
sions. The PAH rovibrational emission is represented by the 
dot-dashed line. Dotted lines show the PAH rotational emis- 
sion. The black model is model A, and the grey model is model 
B. See text for details. 



6. The case of the molecular cloud G1 59.6-1 8.5 

We present here a comparison of our model results with obser- 
vations of the anomalous microwave foreground and IR emis- 
sion. We show how these data can be explained with a coherent 
description of the IR and rotational emission of PAHs, and also 
discuss how the size distribution of PAHs can be constrained. 
The COSMOSOMAS experiment has delivered maps of large 
sky fractions with an angular resolution of 1 degree in 
four freque ncy bands: 10.9, 12 .7, 14.7, and 16.3 GHz. In 



this survey, IWatson et al.l ( 120051) studied the region G159.6- 
18.5 located in the Perseus molecular complex. Figure [TT] 
shows the spectrum of this region from 3 to to 4000 GHz 
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(COSMOSOMAS, WMAP, and DIRBE data). From mid to far- 
IR, the Perseus molecular complex is dominated by a molecular 
ring (G159.6-18.5) surrounding an Hll region. Unfortunately, 
the COSMOSOMAS resolution does not allow these two 
phases to be se parated. The molecular ring is centred on the 
star HD278942 dAndersson et al.ll2000l: liidge et al.ll2006l) . Its 
stellar wind is responsible for the HII region dug into the parent 
molecular cloud. To describe the emission spectrum of G159- 
18.5, we assume it includes an MC component for the ring and 
a WIM component for the Hll region. We discuss below the 
physical para meters of both pha ses. 

According to lRidge et all (l2006h . HD278942 is a BO V star: its 



radiation field is modelled by a blackbody with T^jf ~ 30,000 
K, scaled to get Go ~ 1.6 at the radius of the ring (Go is ob- 
tained from IRAS data). The radius of the Hll region is approx- 
imately 2.75 pc dAndersson et al. 2000l) and implies a WIM 
density «H ~ 1 cm~^ dRidge et al.l2006 ). We obtained the WIM 



Andersson et al. I d2000h 



temperature from CLOUDY: 7 500 K. 
estimated the density and temperature of the gas in the molecu- 
lar ring to be hh ~ 560 cm"^^ and T = (20+ 10) K. For the WIM 
and MC phases, we derived the abundances of electrons, pro- 
tons and C^, and the gas temperature with CLOUDY. We took 
the tot al column density ( MC+WIM) to be A^// = 1.3 x 10"^ 



cm (iWatson et al.ll2005h . From 12 fim IRAS data and our 



PAH emission model, we estimated the total PAH abundance 
to represe nt (240 ± 40) ppm of carbon. 

Following lWatson et al.l ( l2005b . we took the thermal dust emis- 
sion to be a grey body at 19 K with an emissivity index 
y8 = 1.55. Our best-fit parameters for the spinning PAH con- 



tribution (assuming a size distribution 11,(0) ~ a ^■^) are then 

-\21 ^^-2 



cm ^ for 
for the 



(model A) Nc ^ 50 - 216 and Nh = 4.6 x 10 
the MC, Nc ^ 24- 216 and A^^ = 8.1 x 10^' cm 
WIM. In both media, we took m = 0.6 D (Eq.|4]i. Larger PAHs 
are required in the MC, possibly as a co nsequence of grain- 
grain coagulation. However, the results of iDupac et al. (2003) 
indicate that emission of dust around 20 K requires a higher 
emissivity index, close to 2. When adopting 13-2 and a dust 
temperature of 18 K (model B), we find an equally good fit of 
this region in the far-IR (Fig.fTTTi. In this case, the flux around 
100 GHz is not explained. We speculate here that it may come 
from an additional and less abundant population of small PAHs 
in the MC component with: Nc - 24 - 50, m = 0. 1 D and 
Nh - 1.25 X 10-' cm"^ (for the larger PAHs in the MC phase, 
Nh = 4.25 X 10"' cm"^). The corresponding spinning emis- 
sion is shown in Fig. [TT] (moc/e/ B). The PAH si ze distribution 



would thus be bimodal, as already suggested bv lLe Page et al 



(l2003b . As can be seen in Fig. [TT] both models provide a good 
fit to the data. Discrimiting tests of these scenarios will be soon 
possible with the Planck polarized data. Small grains or PAHs 
are expected to be poorly aligned (Lazarian & Draine 2000, 
Martin 2007), whereas big grains are fairly well-aligned with a 
polarization fraction of 5 to 10%. 

7. Summary 

The Planck and Herschel data will soon reveal the emission 
of interstellar dust at long wavelengths. Thanks to their small 
size, interstellar PAHs spend most of their time at low inter- 



nal energy and can spin at frequencies of a few tens of GHz. 
The emission of PAHs is therefore expected to make a signifi- 
cant contribution at long wavelength. Recent observations have 
shown the existence of a 10-100 GHz emission component (the 
anomalous foreground), related to the smallest dust grains. As 
suggested by DL98, the anomalous foreground may trace the 
emission of spinning PAHs. 

In this work, we built the first model that coherently describes 
the emission of interstellar PAHs from the near-IR to the cen- 
timetric range and focused on the long wavelength part of this 
emission. To do so, we derived the internal energy of isolated 
PAHs down to low energies. We included low-frequency vibra- 
tional bands (A > 20 jum), through which PAHs cool at inter- 
mediate to low energies. They are important for their rotational 
emission. In the cooling cascade that follows the absorption 
of a stellar photon by a PAH, we treated rovibrational transi- 
tions in a quantum approach and examined the possibility that 
the hypothesis of vibrational redistribution (IVR) is not always 
fulfilled. Purely quantum effects (recoil due to photon emis- 
sion; transitions that do not change the angular momentum or 
g-bands) are thus naturally included in our description of the 
rovibrational cascade. We obtained the rotational emission of 
PAHs from a random-walk formalism, including all processes 
participating in excitation and damping, namely, rovibrational 
and pure rotational transitions and gas-grain interactions. 
We have shown that the rovibrational emission of PAHs above 
~ 3 mm does not depend on the intensity of the radiation field 
(represented by Go), unlike the mid-IR part of the spectrum that 
scales linearly with it. In the diffuse interstellar medium, PAHs 
may contribute up to 10% of the dust emission around 100 
GHz. We also found the rotational emissivity of PAHs is dom- 
inated by small species (bearing less than 100 C atoms) and 
is hardly sensitive to Go over the range 0.1-100. Using plau- 
sible PAH properties, our model can explain both the IR and 
microwave emissions of a molecular cloud in the Perseus arm 
where the anomalous foreground is conspicuous. The level, 
peak position and width of the rotational emission spectrum 
depend on the fraction of small PAHs, the dipole moment dis- 
tribution (m-factor), and the physical parameters of the gas 
phases present along the line of sight. A quantitative descrip- 
tion of the emission of spinning PAHs will therefore involve 
observations at IR and radio wavelengths (21 cm and contin- 
uum). Comparing the rotational excitation rate obtained from 
our quantum treatment of the rovibrational cascade to former 
works, we showed that the classical approximation used so 
far has little effect on the rotational emissivity (the peak value 
varies by at most 15%). Similarly, departures from the IVR hy- 
pothesis lead to similar emissivity changes. We therefore con- 
clude that a classical description of rovibrational transitions 
and the IVR hypothesis are good approximations for describing 
the rotational emission of PAHs. 

Our results on the influence of the radiation field intensity led 
to a specific prediction that can be tested observationally. If the 
anomalous microwave foreground comes from spinning PAHs, 
it is expected to be correlated with the dust emission in the 12 
jum-IRAS band, mostly carried by PAHs. In regions where Go 
varies significantly, this correlation should improve when the 
12 jum flux is divided by Go, indeed, the IR emission of PAHs 
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scales with Go, whereas their rotational emission is indepen- 
dent of Go- This prediction was tested in a co mpanion paper 
with WMAP and IRAS data dVsard et al.1120091) . 
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Appendix A: Mid-IR absorption cross-section 



In Table lA.n we give the parameters defining the mid-IR vibrational 
bands considered in this work. In all cases we assume their profile to 
have a Drude shape. The bands at 3.3, 6.2,7.7, 8.6, 11.3, and 12.7 yum 
are defined as follows. For the PAH catio ns, we start from the inte- 
grated cross-sections, crAv, of IPech et al. I |2002i) that have been de- 
rived from laboratory data. The corresponding band profiles, however, 
do not provide a detailed match of observations. We therefore use 
band positions and widths as deduced from fits of ISO-SW S spec- 
tra of a number of interstellar regions IVer straete et al J ^OOll ). As in- 
dicated by these observations and others dPeeters et al.ll2002l) . we in- 
clude a broad band at 6.9 yum and split the 7.7 yum into three sub- 
bands at 7.5, 7.6, and 7.8 /im, where we use the observed crAv of 
each sub-band as weights in defining their integrated cross-sections. 
In addition, we introduce a band at 8.3 yum to fill the gap between 
the 7.8 and 8.6 yum bands, and multiply the 8.6 yum band by a fac- 
tor 3 to match observations. F or neutral PAHs, w e use the labora- 
tory integrated cross-section of I Joblin et al. I i ll995l) . and assume the 
same band profiles as for the cations. Furthermore, spectroscopic data 
(ISO-SWS, Spitzer -IRS and UKIRT) reveal other bands at 5.25 and 
5.75 yum, which have been ascribed to combinations of PAH vibra- 
tional modes involving the 1 1 .3 yum band and IR-forbidden modes a t 
9.8 and 1 1.7 /jm, respectively l lRoche et al.ll996l ; lTriDathi et al.l200ll) . 
For the 5.25 yu m, we use the width and intensity ratio to the 11.3 yum 
band given in iRoche et ai] ( Il996l) . The 5.75 yum ba nd has been de- 
rived from the observed spectrum of the Orion Bar jVerstraete et alj 
I2OOII) . We also add the 17.1 fxm band recently seen in Spitzer data 
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Table A.l. Mid-IR bands of interstellar PAHs adopted in this 
work for cations and neutrals. 





V; 


Avi 


CTilNn 


a-i/Nn 


Type^ 


(pm) 


(cm ') 


(cm ') 


(10 ^"cm^) 


(10 -"cm-) 










cations 


neutrals 




3.3 


3040 


39 


2.44 


10.8 


ip 


5.2 


1905 


23 


0.58 


0.58 


op 


5.7 


1754 


60 


0.49 


0.49 


op 


8.3 


1205 


63 


1.74 


1.74 


ip 


8.6 


1162 


47 


5.34 


0.51 


ip 


11.3 


890 


18 


17.3 


18.3 


op 


ion 


ojU 


iU 


J.V 1 


i.l / 


op 


12.7 


785 


16 


5.06 


4.06 


op 


A, 




Av, 


o-.INc 


(Ti/Nc 


Type 


(jim) 


(cm-') 


(cm-i) 


( 10-20 cm^) 


(lO-'^cm^) 










cations 


neutrals 




6.2 


1609 


44 


2.48 


0.52 


ip 


6.9 


1450 


300 


0.40 


0.40 


ip 


7.5 


1328 


70 


2.70 


0.12 


ip 


7.6 


1315 


25 


1.38 


0.06 


ip 


7.8 


1275 


70 


2.70 


0.12 


ip 


16.4 


609 


6 


1.83 


1.83 


ip 


17.1 


585 


17 


2.48 


2.48 


ip 



In-plane (ip) or out-of-plane (op) character of the bands (see Sect. 



dSmith et alj2004 [Werner et al .'2004') an d rec ognized as arising from 
PAHs jPeeters et al.lf2004b ; .Smith et al.„2007l) . 




E (cm-') 

Fig.B.l. Cumulative distribution of vibrational modes for 
coronene (C24//12 ), ovalene (C32//14 ) and circumovalene 



iC(,bH2o), from the iMalloci et al.l (120071) database (solid lines) 



and with our Debye model (dotted lines). 



Appendix B: Vibrational modes and density of 
states of interstellar PAHs 

B.1. Vibrational modes 

Symmetric top (Dg;, symmetry) type PAHs, with Nc carbon atoms and 
Nu = s/bNc hydrogen atoms, have 3(Nc + Nh - 2) vibrational modes 
that can be divided into the following types: (Nc - 2) out-of-plane 
(op) CC modes, 2(Nc - 2) in-plane (ip) CC modes, out-of-plane 
CH bending modes, Nh in-pla ne CH bending mo des, and Nh CH 
stretching (st) modes. FolIowing lDraine & il Eooih . we approximate 
the mode spectrum of each type of vibration with a two-dimensional 
Debye model of maximum energy k@, where is the Debye temper- 
ature. We derive the mode spectra from the following expressions: 



- for the CC modes: 



h(ji>i = k@i 



N, 



for ;■ = l,N, 



(B.l) 



where t = op, ip is the type of mode and A', is the number of CC 
modes of a given type with @„p = 876 K and 0,p = 2318 K with: 

S°" = 3/2 for = 2, 3 and 6'." ■ 
= 3/4 otherwise 



1 for = 2, 3 
1 /2 otherwise 



- for the CH modes: 



ficOj 



IN, 



(B.2) 



(B.3) 



where t = op, ip, st, and A', is the number of CH modes of a given 
type with 0„p = 1281 K, 0;^ = 1672 K and 0„ = 4375 K. 



To this mode spectrum, we add the mode at lowest energy as described 
in Sect. l2.1l Figure lBTTjshows that our Deb ye mode spectra agree with 
that derived from the lMalloci et al. U2007h database. 



B.2. Density of states and degeneracies 

We estimate the PAH density of harmonic vibrational states from 
direct counts, a mode-to -mode convolution method proposed by 
iBever & Swinehar^ ( Il973l) . We start from the rotational density of 
states, classically given by Prot (E) = Vlf/B^^^ with E and B in 
cm- and for a symmetric top ( lBaer&Has3ll996h . We then obtain 
the rovibrational density of states p(E) with the convolution method. 
The molecule's zero-point energy has been chosen as the zero of the 
energy scale, and the calculation was made for bins with finite width 
of 1 cm-'. To calculate the internal energy distribution of PAHs P(E), 
we grouped this very large number of points into broader energy bins 
[E""" , E"""'] with / = 1 to 500. Each energy bin / thus contains many 
states and its degeneracy, gj, is estimated as (with gi = 1) 



p(E)dE (i > 1) . 



(B.4) 



Appendix C: Lines intensities in the rotational 
bands 

In the case of symmetric top molecules, we give below the Honl- 
London factors for the angular part of rovibratio nal transition rates , 
as well as the corresponding transition energies ( lHerzberdll968abh . 
The frequency of the vibrational mode i is noted v/o. 
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- for AJ = 0,-- 
transitions): 



1 and AK = (parallel transitions or out of plane 



for AJ 
tions): 



2{(J+lf-K^) 
(J + l)(2J + 1)2 
2(J^ - K^) 
J {2 J + 1)2 
2K^ 

J(J+l)(2J+ 1) 



0,±1 and AA: 



and 



and 



and 



v,o - 2B(J + 1) 



y,o + 2BJ 



(C.l) 

:1 (perpendicular or in plane transi- 



_(J + 2± K)(J +\±K) 

~ 2(y + i)(2i+ 1)2 

= v,o - 2B{J + 1) + (B - C)(l ± 2K) 

^ _ (7 - 1 + K){ J + K) 
~ 27(27 + 1)2 

vT = y,o + 267 + (S - C)(l ± 2K) 
(7 + 1 ± K){J + K) 



^KJO 



27(7 + 1)(27 + 1) 



v,o + {B-C)(l±2K) 



Appendix D: Rate of angular momentum change 
for gas-grain interactions 

We describe below the gas-grain interactions considered in our model 
of PAH rotation and the rate of change of 7 they induce, (t"'A7). 

D.1. Collisions with gas atoms and plasma drag 

For the collisions with gas neutrals and ions, as well as the plasma 
drag, we apply the results of DL98 to planar PAHs. In the case where 
7 >> 1, we use the correspondence principle to write hj = 1^.0) and 
obtain the following rates: 



(t-'A7) = -5.2 X 10" 



V lOOKUOOcm-3/ 



xJxF 



(D.l) 



for the damping contribution, and 



(r-'A7) = 3.9 X 10^ (^)' - ] - ^ (D.2) 

for the exciting contribution, where F and G are normalized rates 
defined in the Append ix B of DL98. We use the formalism of 
iBakes & Tieleii^ j 19941) to estimate the average charge of PAHs of 
a given size. 



that Ht molecules are ejected from the edge of the PAH with a cosine 
law: 



(T-'A7)„, 



3.4X10- (^f ( ,) 
V 50 / V 100 cm-3 / 



100 K \ 1.5 eV 



(1-/h) s- 



(D.3) 



The above numerical values are based on the following assumptions: 

(i) all molecular hydrogen is formed on PAHs with [C/H]pah = 

4x 1 0~^ a Hj formation rate Rf = 3 x 10-"(r/70K)'/2 cm^s"' jjurj 
1 19751) and E^j = 1.5 eV ; (ii) the distribution of formation sites has 
an asymmetry of 1 site and we assume that the site in excess is al- 
ways at the same location on the molecule ; (iii) we neglect the in- 
fluence of cross-over even ts that may reduce the angular momentum 
jLazarian & DrainellT999h . With all these assumptions, the spin rate 
due to H2 formation estimated here is an upper limit. This rate, how- 
ever, remains small compared to the other gas-grains processes. 



D.3. Photoelectric effect 

Stellar UV photons can pull out electrons from grains. These photo- 
electrons carry away a significant kinetic energy (~ leV) that heats 
the interstellar gas and impulses grain rotation. If we assume that the 
photoelectrons leave the grain surface with a cosine law distribution, 
we have 



(T-'A7)„ 



t;; xO.15 



(D.4) 



where t^I is the rate of photoelectro ns ejections calculated with the 
formalism of lBakes & TielensI i ll 9941) {tJ,1 = 1.6 x 10"** s"' for Nc = 
50 in the DIM). 



Appendix E: Rotational temperature 

A rotational temperature Tj can be defined from the Maxwell distri- 
bution of angular momentum, «(7) = no exp(-72/7Q): 



r,(K) : 



BJl 



(E.l) 



Figure lETU shows this temperature as a function of PAH size for sev- 
eral interstellar environments. We see that is subthermal in the case 
of the DIM, and suprathermal for MCs and the Orio n Bar. These re- 
sults are in good agreement with lRouan et al. I l ll997h . The rotational 
temperature is a relevant parameter in the study of the width of PAH 
vibrational or electronic transitions. The latter have been proposed as 
the origin of some unidentified diffuse interstellar bands (DIBs). 



D.2. Rocket effect 

Ejection of H or H2 from the edges of PAHs may yield a significant 
rotational excitation if it occurs as ymmetrically, thus g enerating a sys- 
tematic torque by the rocket effect jRouan et al J 19921) . We note below 
Egj the kinetic energy of the ejected fragment. In the case of H2, this 
will happen if this molecule forms on preferential sites by chemisorp- 
tion, and if the distribution of these sites on the PAH is asymmetric as 
a result of dehydrogenation. We calculate the change of 7 assuming 
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Fig.E.l. Rotational temperature for PAH cations as a function 
of their size for several insterstellar environments: CNM (black 
Une), WNM (gray dashed line), WIM (gray line), MC (black 
dotted line) and the Orion Bar (black dashed line). 



